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Abstract 

Bernstein inequalities and inverse theorems are a recent development in the 
theory of radial basis function(RBF) approximation. The purpose of this 
paper is to extend what is known by deriving LP Bernstein inequalities for 
RBF networks on M. d . These inequalities involve bounding a Bessel-potential 
norm of an RBF network by its corresponding LP norm in terms of the 
separation radius associated with the network. The Bernstein inequalities 
will then be used to prove the corresponding inverse theorems. 
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1. Introduction 

When analyzing an approximation procedure, there are typically two es- 
timates that one is interested in determining. The first is a direct theorem 
that gives the rate at which a function may be approximated in terms of its 
smoothness, and the second is an inverse estimate that guarantees a certain 
amount of smoothness of an approximant based on its rate of approximation. 
Both are equally important, and if the results match up appropriately, they 
can be combined to completely characterize smoothness spaces in terms of 
the approximation procedure. In this paper, our goal will be to prove an 
inverse theorem for RBF approximation. The usual way one does this is by 
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first deriving a Bernstein inequality, since a standard technique can then be 
applied to prove the inverse theorem. 

Bernstein inequalities date back to 1912 when S.N. Bernstein proved the 
first inequality of this type for L°° norms of trigonometric polynomials, [l[ . A 
generalization can be found in j^j; this result, which is credited to Zygmund, 
states that any trigonometric polynomial T of degree n satisfies 

||T (r) || <n r \\T\\ 
II I \ P — ii Hp 

for 1 < p < oo. However, the first example of a Bernstein-type inequality for 
RBF approximants was not proved until 2001, Q. Then in 2006, Narcowich, 
Ward, and Wendland derived a more standard type of Bernstein inequality, 
0] . They proved L 2 Bernstein inequalities for approximants coming from an 
RBF approximation space on M. d where the Fourier transform of the RBF 
has algebraic decay. In the same year, Mhaskar proved LP Bernstein inequal- 
ities for certain Gaussian networks on W 1 , jjjj], and lastly, in Q, Mhaskar, 
Narcowich, Prestin, and Ward were able to prove Bernstein inequalities in 
LP norms for a large class of spherical basis functions (SBFs). 

In this paper we will be concerned with RBF approximants to functions in 
L p (M. d ) for 1 < p < oo. The approximants will be finite linear combinations 
of translates of an RBF $, and the translates will come from a countable set 
X C M d . The error of this approximation, which is measured in a Sobolev- 
type norm, depends on both the function $ and the set X. Therefore, given 
an RBF $ and a set X, we define the RBF approximation space Sx(&) by 




By choosing $ and X properly, one is able to prove results about rates of 
approximation as well as the stability of the approximation procedure. 

Our goal will be to establish LP Bernstein inequalities for certain RBF ap- 
proximation spaces Sx($), and these inequalities will take the form | \g\ \ Lk>p < 
Cq x k \ \g\\ p , where L k ' p is a Bessel-potential space. To prove this, we will use 
band-limited approximation with the bandwidth proportional to l/qx- Thus 
l/qx acts similarly to a Nyquist frequency, and viewing l/qx as a frequency, 
we can see the connection to the classical Bernstein inequalities for trigono- 
metric polynomials. In particular, bandwidth is playing the role of the degree 
of the polynomial from the classical inequality. 
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The RBFs that we will be concerned with have finite smoothness; exam- 
ples include the Sobolev splines and thin-plate splines. The sets X will be 
discrete subsets of M. d with no accumulation points. For the inverse theorem, 
we will additionally require that there do not exist arbitrarily large regions 
with no point from X. 

1.1. Strategy 

The basic strategy that we will use is the following, which is the same 
as the one used in [4J. Given g = ^^ gl a^$(- - () 6 Sx($), we choose an 
appropriate band-limited approximant g a , and we have 

||fllli>,p — ||<?o-||i,fc,i< + \\g — 9a\\L k 'P ■ 
We then split the second term into two ratios. 

ii ii ^ ii ii (W a \\p\\9 ~ 9a\\ L k, P \ M 1 1 m 

\\g\\ L k,p < \\9a\\ Lk , P + jTT: rr-r: \\g\\ p . (1) 

\\\9\\p W a Wp J 

The term ||a|| / ||<7|| p will be bounded by a stability ratio 7Zs,p that is in- 
dependent of the function g. We will then need to bound the error of ap- 
proximating g by band-limited functions. Combining these results with a 
Bernstein inequality for band-limited functions, we will be able to prove the 
Bernstein inequality for all functions in S$(X), and afterward the corre- 
sponding inverse theorems will follow easily. 

1.2. Notation and Formulae 

For any approximation procedure, one would like to determine the error of 
the approximation and the stability of the procedure. When considering an 
RBF approximation space Sx($), these quantities are bounded in terms of 
certain measurements of the set X. The error of approximation will typically 
be given in terms of the fill distance 

hx = sup inf \\x — £|| 2 , 

which measures how far a point in M d can be from X, and the stability of 
the approximation will be determined by the separation radius 
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which measures how close two points in X may be. In order to balance the 
rate of approximation with the stability of the procedure, approximation will 
typically be restricted to sets X for which hx is comparable to qx, and sets 
for which the mesh ratio px '■= hx/qx is bounded by a constant will be 
called quasi-uniform. 

Many of the results in this paper will be proved by working in the Fourier 
domain, and we shall use the following form of the Fourier transform in R d . 



The function spaces that we will be mainly interested in are the Bessel- 
potential spaces L k ' p (M. d ), which coincide with the standard Sobolev spaces 
W k ' p (M. d ) when A; is a positive integer and 1 < p < oo, cf. jsl, Section 5.3]. 
The Bessel potential spaces are defined by 

L k, P = {f: f ={1 + \\.\Q-»/% g e L p (R d )} 
for 1 < p < oo, and they are equipped with the norm 

li*.p = WdWp ■ 



1.3. Radial Basis Functions 

In order to prove the Bernstein inequalities, we will need to require 
certain properties of the RBFs involved. As much of the work will be 
done in the Fourier domain, we state the constraints in terms of the RBFs' 
Fourier transforms. Given a radial function $ : M. d — > R with (general- 
ized) Fourier transform $, let : (0, oo) — > R be the function defined by 
<b(oj) = (f)(\\u\\ 2 ). We will say a function $ is admissible of order /3 if there 
exist constants C\, C 2 > and > d such that for all a > 1, x > 1/2, and 
l<k-= [(d + 3)/2], we have 

(i) Ci< (f>{(Jx){l + {axff 2 <C 2 



((l + (ax)y/ 2 0(ax)) (O 



< Co 



Two particular classes of admissible functions are the Sobolev splines and 
the thin-plate splines. The Sobolev spline $ of order > d is given by 

® = C\\-\\t d)/2 K {d ^y 2 (M 2 ), 
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where K is a modified Bessel function of the third kind. This function pos- 
sesses the Fourier transform 

$ = (l + ||.||5)-# 2 . 

For a positive integer m > d/2, the thin-plate splines of order 2m take the 
form 



1 2m— d r ij 

| 2 a odd 

| 2 log ||-|| 2 a even 



and possess the generalized Fourier transforms 

$ = c||-||- 2m . 

1.^. Bessel Functions and Fourier Integrals 

A cZ-dimensional Fourier integral of a radial function reduces to a one 
dimensional integral involving a Bessel function of the first kind. In what 
follows, we will discuss some of the properties of these Bessel functions as 
well as some propositions which will be useful for bounding Fourier integrals. 



Proposition 1.1. ([10L Proposition 5.4 & Proposition 5.6]) 



(1) £{Z»UZ)} = ?Ju-l{z) 



(2) J 1/2 (z) = J^sm(z), J- 1/2 (z) = J^cos(z) 



(3) J u (r) = yj£ cos(r - if - f ) + 0(r- 3 / 2 ) for r -»■ oo and v G R 

(4) Jf /2 {r) < ^ forr > and I G N 

^ hm r ^o r- l Jf /2 {r) = 2 i v ^ l/2+l) for I en 

The next proposition makes use of integration by parts in order to bound 
the Fourier integral of a function whose support lies outside of a neighborhood 
of the origin. 

Proposition 1.2. Let a > 1, and let f G C n ([0, oo)) for some natural 
number n > 1. Also, assume there are constants C, e > such that f = on 
0, |] and < Ct~ d ~ e for j < n. Then there is a constant C e such that 

Jy 2 f(t)t d / 2 J (d _ 2)/2 (at)dt <C € a~ n 
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Proof. We first define a sequence of functions arising when integrating by 
parts. Let f = f,f 1 = /', and = for j > 2. Note that when j > 2, 

there are constants cjj such that fj(t) = Yll=i c j,if^ l Kt)t^ 2j+l+1 , and therefore 
\fj(t)\ < Ct~ d ~ j+1 ~ e . Applying the Dominated Convergence Theorem, we 
have 



1/2 



Ht)t dl2 J {d -2)/2{cd)dt 



lim 

b— ¥OD 



1/2 



h{t)t d/2 J {d - 2 )/2{at)dt 



After integrating by parts and taking the limit we get 

hit) 



1/2 



fo(t)t d/2 J id ^)/2(at)dt 



1 

a 



1/2 



t 



t d/2+1 J d/2 (at)dt 



Integrating by parts j times, we have 



fo{t)t d l 2 J [d _ 2)/2 {at)dt 



1/2 



and therefore 



= — r lim 



" fj(t) td/ 2 +jj l{at)dt 
1/2 ■ 



t 



f j+ i(t) 



1/2 



t d/2 +j+ lj d/2+ .( at)dt 



POD 

/ h(t)t d l 2 J (d _ 2)/2 (at)dt 

Jl/2 





POD 








Jl/2 



fn(t) 



t d/2+n J d /2 + n-Mt) 



dt. 



□ 



We will also be need to bound Fourier integrals of functions which are 
identically one in a neighborhood of the origin and have compact support. 
The next proposition addresses this. As in the previous case the proof relies 
on integration by parts. 

Proposition 1.3. Let a > 1, and let f be a function in C n ([0, oo)) for some 
natural number n > 1. Also, assume f = 1 in a neighborhood ofO and f = 
for x > 2. Then there is a constant C such that 



J,tf(t)t d / 2 .J {d ^ 2)/2 (at)dt 



< Ca~ n 
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Proof. We first define a sequence of functions arising when integrating by 

parts. Let f (t) = tf(t), fx = and f 3 - = for j > 2. Note that 

when j > 2, fj(t) = 0(t~ 2 i +2 ) as t — > 0. After integrating by parts n times, 
we have 



fo(t)t d/2 J { 



(d-2)/2 



at)dt 



1 






I 



^t d / 2+n J d/2+n ^(at)dt 



□ 



1.5. Sums 

In the proof of the Bernstein inequalities, we will show that the constants 
involved are independent of the number of centers. The only dependence of 
the constants on the set X will be in terms of the separation radius qx- A 
key result for showing this independence is the following proposition. It deals 
with bounding a sum of function values on a finite set. 

Proposition 1.4. Let X C M d be a countable set with qx > 0, and let 

Y = {yj}f =1 be a subset of X such that \\yj\\ 2 > qx f or 1 < j < N. If 



f- 



— > 



is a function with \f(x)\ < C \\x\ 



for some C, e > 0, then 



N 



Ei/(%-)i<3 rf (i 



l/e){Cq x d ^) 



Proof. We can bound the sum using the volume argument found in the proof 
of 10, Theorem 12.3]. Following the same procedure, we define 

E m = {x E R d : mq Y < \\x\\ 2 < (m + l)qy}- 

for each positive integer m. Now by comparing the volume of E m to the 
volume of a ball of radius qy, one finds that j^Y D E m < 3 d m d ~ 1 . Therefore, 



N 



A' 



c 



d+e 
3 =1 nojn2 



!\\vi\\\ 



< 



< 



V(#^nE m ) max 

m=l 



C 



\y\\2 



oo 



c 



< 3 d (l + l/e)Cg 



-d- 



□ 
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2. Stability 



One of the essential results for proving the Bernstein inequalities is a 
bound of a stability ratio for S x (&). We define the L p stability ratio 1Zs, p 
associated with this collection by 



TZ s ,p = sup 



a 



p — qujj ! ! ^ , 

lp 



Sx(*)9S^0 | IP | 



where (7 = Xlgex a £^K" ~~ 0- The goal of this section is to bound the stability 

ratio by Cq d Jr V 13 for some C independent of a and X, where p' is the conju- 
gate exponent to p. For this section, we will assume we are working with a 
fixed countable set X C M. d with < qx < 1 and an admissible function <3> 
of order f3. 

To begin, fix Y = {£j}f =1 C X and 5 = J^li a j^(' — £?')■ We wm derive 
a bound for ||a|L/ ||<?|L and show the bound is independent of Y and a. 
The strategy for proving this is as follows. Let K be a smooth function and 
define K a (u>) = K(ou/a). We will then consider the convolutions K a *g(x) = 
J2j f =i a j^(r * ®(x — For an appropriate choice of a, the interpolation 
matrix (Ar)j,j — K a * — £j) will be invertible, and the norm of its inverse 
will be bounded. Then a = A~ l (K a * g)\y and ||a|| < | j^4~ x 1 1 \\K a * g\y\\ p - 
We will then be left with bounding \\K a * g\y\\ p in terms of qx and \ \g\\ p . 

2.1. Convolution Kernels 

We now define the class of smooth functions with which we will convolve 
g. Let /Ci be the collection of Schwartz class functions K : R d — > H. that 
satisfy: 

(i) There is a k : [0, 00) — > [0, 00) such that K(u) = k(||o;|| 2 ). 

(ii) n(r) = for r e [0, 1] and k is non-vanishing on an open set. 

For each K e /Ci, we define the family of functions {K a } a >i by K a (uj) = 
K(u/a). Note that property (ii) requires each function K a to have a Fourier 
transform which is in a neighborhood of the origin, and as a increases, so 
does this neighborhood. The convolution $ * K a will retain this property 
and allow us to obtain diagonal dominance in A a . 
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Before moving on, we will need to determine certain bounds on the func- 
tions in JC\. First we need an L°° bound. 



\K a (x)\ < C d / K a (u)du 

/oo 
K(t)t d - l dt, 



so 



\KJx)\ <C K4 a c 



(2) 



Next we will need a bound on K a for r = \\x\\ 2 > 0. Writing K a as a 
Fourier integral, we see 



\K a (x)\ = r-^' 2 
and by a change of variables, we have 

\ Ka (x)\ = a d/2 + l r -(d-2)/2 

Therefore by Proposition 11.21 



K 



{t/a)t d,2 J {d ^ 2) , 2 {rt)dt 



K(t)t d ' 2 J {d „ 2)/2 (art)dt 



\KJx)\ < C 



a 



d/2+l-l d 



KA 



r (d-2)/2+l d ■ 



(3) 



Using 02]) and ([3]) , we will prove a bound of L q norms of linear combina- 
tions of translates of K a . 

Proposition 2.1. Let T : R N — >■ L 1 (M d ) H L°°(R d ) be the linear operator 
defined by 

N 



T{h) = Y J h j K (T {x-i j ). 

3=1 



Then 



\\T{h)\\ <c K y^' 1 + 



1 



(d-2)/2+Z d 



l/ll 
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Proof. After proving the bound in the cases q = 1 and q = oo, the result 
will then follow by the Riesz-Thorin theorem (cf. P, Chapter 5]). First, by 
Proposition II .41 we have, 



N 



j =1 \\x-£j\\>q 



< C K4 [o d + 



r d/2+l-l d 
(d~2)/2+l a 



Simplifying the previous expression, we obtain 

N 



{ 

J2\Ka(x-Zj)\ <C K>d a d l 1 

.7 = 1 ^ 



( ffq yd-2)/2+l d 



Therefore 



\\T(h)\l 



N 



< 



AT 



< C M a d I 1 4 



oo 

1 



( aq yd-2)/2+l a 



Now in the case q = 1, 



TV 

to iii = 

N 

< CdWIi 



(4) 



□ 
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2.2. Interpolation matrices 

The interpolation matrices (Ar)i,j = K c * — £,) will be shown to be 
invertible by the following lemma. In addition, we will at the same time find 
a bound for the l v norm of A~ l . Given an n x n matrix A with diagonal 
part D and F = A — D we have the following. 

Lemma 2.2. ([4, Lemma 5.2]) If D is invertible and \ \D~ 1 F\\ 1 < I, then 
A is invertible and 1 < H-D -1 ^ (1 - | \D^F\ 

The diagonal entries of A a are equal to K a * $(0), and the off diagonal 
absolute column sums are of the form \K-a * ~ I n order to 

apply the lemma, we must bound the former from below and the latter from 
above. First we have 



Next, we need a bound on \K a * $(x)| for i^O. Since K a * & has a radial 
Fourier transform in L^IR^), we can write it as a one dimensional integral. 
Note that in the following integral r = \ \x\\ 2 . 




and hence 



K a * $(0) > C w & 



(5) 





{o-tf<P{o-t)t d/2 J {d „ 2 )/2{o-rt)dt , 



and by Proposition 11.21 



a 



(d+2)/2-/3 



K a *®(x)\<C*, K>d 



r {d-2)/2( ar y d 
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With this estimate we can bound the off diagonal absolute column sums of 
A a . Using Proposition 11.41 we have 



£ \k„ * m - o)i < c,, K , d {(7qx){d _ 2)/2+ld . (6) 



We are now ready to apply the lemma. Define 

'2C? a V((rf-2)/2+/ d ) 



M = max < 1 



r 1 



where the constants C\ Kd and C\ K d are from ([5]) and ([6]) respectively. We 
then define <7o = M/qx, so 

Therefore 

H^^lli ^ C$ :K;d a d , 

and in terms of qx, 

IKII^C^g^. (7) 

As A ao is self-adjoint the same bound holds for H^^Hoo- The Riesz- 
Thorin interpolation theorem can then be applied to get 

\\4£\\ p <c 9 #j£- p (8) 



for 1 < p < oo. 



2.3. Marcinkiewicz-Zygmund type inequality 

To finish the bound of the stability ratio we require a bound of a discrete 
norm by a continuous one. To accomplish this, we can use an argument 
similar to the proof of [5|, Theorem 1]. Let fi denote Lebesgue measure on 
R d and let v = £f =1 <%.. 

Proposition 2.3. If f E L^D and K e K, 1} then 



\K* * fWWp < C K4 q x 



d/ P 

i i ; i i 
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Proof. First, define T(r, /) = f Rd K m) (- — y)f(y)dr y . Then choose a com- 
pactly supported simple function H such that ll-H" IL , = 1 and 

P i 



N 



\K ao *f\ Y \\ p = j2m)K« *f(& 

3=1 



Appealing to Fubini's theorem, we get 



\\T^J)\\ hu = I T( f i,f)(x)H(x)du :i 

K ao (x - y)f(y)H(x)dfiydu x 
K ao (x - y)f(y)H(x)du x dfj, y 
T(u,H)(y)f(y)d^ 

Now we may apply Holders's inequality and Proposition 12.11 as follows: 

WTinJ)^ < \\T(v,H)\\ p> \\f\\ Pjfl 

N 



3=1 



P ,M 



< Ck4<1x 



-d/ P 



□ 



2-4- Stability Ratio Bound 

We are now in a position to prove the bound on the stability ratio for 
p G [1, oo]. Recall X is a countable subset of M. d with < qx < 1, and <3> is 
an admissible function of order /3. 

Theorem 2.4. Let 7Zs, P be the stability ratio associated with Sx{&)- Then 



K 



S,p 



sup -— 

s x m 39^0 \ \g 



P < r n 1 



d/p'-fS 
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Proof. It has been shown that the interpolation matrix (A aa )ij = K ao *^>(C,i — 
£j) is invertible. Therefore ||a|| p < 1 j^, 1 ) | \\K ao * <7|y|| p . Using (jSJ), we get 

< C^ Ktd q^ \\K ao * g\ Y \\ p - 

Finally, applying the M-Z inequality and taking the infimum over K G K\ 
gives the result. □ 



3. Band-Limited Approximation 

As in the previous section, X will be a fixed countable set with < q x < 
1, and $ will be an admissible function of order (3. At this point, we are left 
with bounding the two remaining terms of ([1]). This will require choosing 
band-limited functions that approximate the elements of <Sx(3>) and satisfy 
the Bernstein inequality as well. In particular, given g G Sx(&) we need to 
find a band-limited function g a so that 



lis - 


■9a 






\a\ 


p 



for 1 < p < oo. Since most of the work will be done in the Fourier domain, 
we will impose the condition k < j3 — d. 

3. 1 . Band-Limited Approximants 

We begin by defining a class of band-limited functions. A function g G 
Sx{*&) will be convolved with one of these functions in order to define its 
band-limited approximant. Let JC2 be the collection of Schwartz class func- 
tions K : ~R d — > R that satisfy the following properties: 

(i) There is a non-increasing k : [0, 00) — > [0, 00) such that K(u) = 

K (IMI 2 ) 

(ii) k(u) — 1 for u < |, and k(u>) = for uj > 1 

For each K G K.2, we define the family of functions {K a } a >i by K a {oj) = 
K(uj/a). Band-limited approximants to g G Sx(&) are then defined by 
g a = K a * g. The first thing we must check is that g a satisfies the Bernstein 
inequality. The following lemma addresses this issue. 



14 



Lemma 3.1. Let f G L 1 (R d ) n L°°(R d ) ; and let K e 1C 2 . Then 

for 1 < p < oo and any positive integer m. 

Proof. We begin with an application of Young's inequality. 



\\f*K a \ 



{i + \\-\\l) m ' 2 kJ 



(i + ||-||^) 1/2 ^(i + ||-||^) (w - 1)/2 ^/"" 



< 



(i 
(i 



^ 1/2 k 2 



(i + \\-\\lY m -^ 2 kJ 

11/ * ^<x|ll / m-l,P • 



•\\2) XV 2<r 

Now it is known that there exist finite measures v and A such that 

(1 + Hxll 2 ,) 172 = v{x) +27r||x|| 2 A(x), 
cf. (§, Chapter 5]. We therefore have 



(1 + I|-Il^) 1/2 ^ 



vK 2a + 27T ||-|| 2 Ai^2<T 



< ||i/*jr 2ff || 1 + 27r 



A * 



I * I la -^2a 



and it remains to prove 



< C 1 



1 2 "^2(7 



I * Ma -^2 CT 



< Ccx. First 



I2 ^2 CT 



< C d \\u)\\ 2 K 2a {u)dw 



2(7 



C d / [ — I t^dt 



C d a l+d I K I - I t fl dt 
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and hence 



Now, for ||x|| 2 = t > 



1 v 



1 2 "^2a 



1 2 -^2 CT 



-(d-2)/2 



l+d 



(9) 



2a 



tn 



A t d ' 2 j { 



2a J 



(d-2)/2 



[rt)dt 



r -(d-2)/2 a d/2+2 



f tK f0 t d / 2 J {d _ 2)/2 {art)dt 



so by Proposition 11.21 



HI2 ^2<r 



< c 



-(d-2)/2 a d/2+2 



(ar) ld 



(10) 



Utilizing inequalities [9] and HUJ we get 



< C d a 



00 ■ / ll»ll 3 >^ 

00 r -(d-2)/2 (T d/2+2 



l 2 ^2a 



dx 



r d - l dr 



err 



□ 



Now if we define <j\ : = 1/qx, then the next corollary follows easily. 
Corollary 3.2. For all g G S x {$), 

\\g * K ai \\ LklP < c<s> tK ,d<ix k WdWp 

3.2. Approximation Analysis 

Now that we know the band-limited approximants to the elements of 
Sx{&) satisfy the Bernstein inequality, we must bound the error of approx- 
imation in L k,p . We will begin by bounding the approximation error in L k)l 
and L k '°° and then use interpolation to obtain the result for all other values 
of p. In both extremal cases, this reduces to bounding the error of approxi- 
mating the RBF by band-limited functions. For p — 1 this is straightforward; 
however, the case p = 00 is more involved. 
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In order to simplify some expressions, we define the functions 



E, 



§,K,k 



((1 + 



|2xfc/2 



($-$* J ftT CT1 ) A ) N 



h(t) := 0(t)(l + t 2 f /2 . 

If we are to bound the error of approximating $ by band-limited functions, 
we will certainly need a point-wise bound of E$ t K,k- Let us begin with an 
L°° bound. 

E*,kjo(x) < [ (l-K ai (io))H^ + \Ml) k/2 doo 

JR d 

(1 - K(t))/l((T 1 t)t d - 1 



d-(3+k 



Therefore 

Next, for ||x|| 2 = r > 



l 1/2 (l/ ffl 2 + *2)(^-fc)/2 



dt. 



(11) 



E$K,k{x) 



-{d-2)/2 



{1 + t 2)(P-k)/2 1 



rt)dt 



<7l 



d/2+l-/3+k 



r (d-2)/2 



(1 - K{t))h{aib) 

1/2 (iM 2 + t 2 )^)/ 2 



t d/2 J, d ^2)/2{(yirt)dt 



Therefore by Proposition 11.21 



£7. 



s (a0 < a 



-(d-2)/2 (7 d/2+l-/8+fc 
((7ir) id 



(12) 



With these results we are now able to bound the error in approximating 
S x {$)- Let y = be a finite subset of X, and let K E JC 2 



Theorem 3.3. Given g = J2f=i a j^(' ~ £j) e Sx(&), we have 

II ts 1 1 ^ /^» /3—k—d/p' II || 

||#-#* <C 9>d q x \\a\\ 



for 1 < p < oo. 
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Proof. We will show this holds when p = 1 and p = oo, and the result will 
follow from the Riesz-Thorin interpolation theorem. Letting g ai = g * K ai , 
we have 



\9 floi llz,fc,°° 



N 



N 



J^k,oo 



3=1 3=1 
N 

3=1 



and using the notation for E§ % & above 



N 



I \9 ~ 9cti 1 1 L fc .°° — I l a l I 
Now by (ITT]) . ( IT2]) . and Proposition II .41 we have 



N 



3 =1 ||a5-fill 2 >«' 

1 



-/M-fc 



1 



(d-2)/2+Z 



For p = 1, we have ||g — 11^,1 < \\a\\\ ll^^^lli; an d 



l^.x^Hi < Crffe I l-E^K,*;! loo + / 

J\\x\\ 



E. 



<!>,K,k 



x)dx 



2><2X 

oo d/2+l-j8+fc-Id 
r (d— 2)/2+Z d — d+l 



dr 



We finish the proof by taking the infimum over K e /C 2 - D 

4. Bernstein Inequalities and Inverse Theorems 

In approximation theory, there are a variety of applications for Bernstein 
inequalities. While they are most commonly associated with the derivation 
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of inverse theorems, they can also be useful in proving direct theorems. For 
example, a Bernstein inequality for multivariate polynomials is used in cer- 
tain RBF approximation error estimates, cf. [10!, Chapter 11]. However, in 
this paper, we will only address the Bernstein inequalities themselves and 
their matching inverse theorems. 

With the bound of the stability ratio and the band-limited approximation 
estimate in hand, we are in a position to prove the Bernstein inequalities. Let 
X be a countable set with < qx < 1, and let $ be an admissible function 
of order /3. 

Theorem 4.1. If k < (3 — d, 1 < p < oo, and g 6 S x (&), then 

\\g\\ L k, P < C$ td q x k \ \g\\ p 
Proof. Let g ax be the previously defined approximant of g. Then 

IMiik.P — Hfi'crJIife.p + \\g — ^crillife.p 



a 


\p 


\g - 


9<n\ 


\L k -P 


IM 


\p 


a 


\p 



\9n\\Lh*+ TTTi iTTl 110 



Applying Theorem \2A\ Corollary 13. 2^ and Theorem 13. 3[ we get 
II II s II II ( H Q IUlff -9<7i\\l** \ ii v 

\\9\\L k -p — iiy<riii£fc.p "I 7i n n n I \\9\ 



< C^ d q x k \\g\\ p + [C* td q$ ) [C^ d q x * p ' ) ||^| 

< C^ d q x k \\g\ 



v 



□ 



Having established Bernstein inequalities for Sx(&), we can now prove the 
corresponding inverse theorems. In what follows, X n will denote a sequence 
of countable sets in M d such that X n C X n+1 , p Xn = p n < C, and < h n , q n < 
2~ n , and $ will be an admissible RBF of order /3. We will additionally use the 
notation S n = Sx n ($)- In this situation, we define the error of approximation 
by 

£(f,Sn) Lp = inf \\f - g\\ L{Rd) - 

Using the standard technique for proving inverse theorems from Bernstein 
inequalities, we will show that if a function is well approximated by S n then 
it must lie in some Bessel-potential space. 
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Theorem 4.2. Let 1 < p < oo, and let f G L p (R d ). If there is a constant 
Cf > 0, independent of n, and a positive integer I such that 



£(f,S n ) Lp <c f h l n , 
then for every < k < min{/3 — d, I}, f G L k,p . 

Proof. Let f n G S n be a sequence of functions satisfying ||/ — f n \\ p < 2cfh l n . 
Note that f n G S m for m > n. We now have 

\\fn+l — fn\\ L k, P < C^dl — J — /n|L 

< ^,^l(ll/n + l-/|| p +||/-/n|| p ) 

< c^,djK+i(Ki+i + Ki)i 



and since h n < 2 it follows that 

\\fn+l — fn\\ L k, P < C*,d,/2 ^ fc " )n . 

This shows /„ is a Cauchy sequence in L k < p . Since L fe ' p is complete, f n 
converges to some function h G L k ' p . Since /„ converges to both / and h in 
L p , f = h a.e., and therefore / G □ 
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